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Abstract 



We present calculations on the deformation of two- and three-layer electret sys- 
tems. The electrical field is coupled with the stress-strain equations by means of 
'^ ■ the Maxwell stress tensor. In the simulations, two-phase systems are considered, 

and intrinsic relative dielectric permittivity and Young's modulus of the phases are 
CJ I altered. The numerically calculated electro-mechanical activity is compared to an 

analytical expression. Simulations are performed on two- and three-layer systems. 
Various parameters in the model are systematically varied and their influence on the 
resulting piezoelectricity is estimated. In three-layer systems with bipolar charge, 
^^ ' the piezoelectric coefficients exhibit a strong dependence on the elastic moduli of 

■^ ■ the phases. However, with mono-polar charge, there is no significant piezoelectric 

effect. A two-dimensional simulation illustrated that higher piezoelectricity coeffi- 
■r^ij- ' cients can be obtained for non-uniform surface charges and low Poisson's ratio of 

^^ ' phases. Irregular structures considered exhibit low piezoelectric activity compared 

"^ . to two-layer structures. 
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. rr , 1 Introduction 



Charged dielectric films with soft elastic properties and at least one free 
(undamped) surface can be used as electro-mechanical and electro-acoustic 
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Fig. 1. Computation domain for the two-layer system. 

transducersjll, 1^. Traditional materials for such applications have been inor- 
ganic crystalline substances [3j. The electro-mechanical effect (piezoelectricity) 
observed in organic polymeric materials differs in several aspects from that of 
traditional inorganic piezoelectric materials. Apart from morphological differ- 
ences, piezoelectricity in polymers is based either on (i) oriented molecular 
dipoles (domains) or on (ii) trapped charges whose distribution breaks the 
symmetry inside the non-uniform anisotropic material. The underlying mech- 
anisms of the phenomenological piezoelectricity in these materials have pre- 
viously been presented in the literaturej^y, |6|,|7[. Recently, several potential 
candidates for large piezoelectric effects in heterogeneous or porous poljTiier 
systems are reported |sl, Isl llO, lulll^, ll^,ll4- Although there exist analytical 



models for modeling simple geometries and uniform charge distributions, nu- 
merical solutions may be preferable for complex geometries [15. 16 . 17] and non- 
uniform charge distributions. In this paper, the finite-element (FE) method is 
first apphed to solve the electric field in two- and three-layered structures with 
interface, surface or volume charges. We have considered non-uniform charge 
distributions and irregular geometries as two-dimensional cases. The obtained 
electric field is later used in calculating the displacement vector by consider- 
ing the Maxwell stress tensor. The results of the one-dimensional simulations 
illustrate that the use of the Maxwell stress tensor generates a perfect agree- 
ment with the analytical model. Since there are no analytical expressions for 
arbitrary geometries and structures with several materials and space-charge 
regions, the application of the FE method to such problems results in a better 
understanding of composite properties that may lead to new tailored materi- 
als. 



2 Numerical Modeling 



A model for the electro-mechanical response of a double-layer dielectric system 
with an interface charge p at the double-layer boundary has been presented 
by Kacprzyk et al. p]. Layered systems have been shown as the optimal matrix 
microstructures for piezocompositesflSl]. In the model, the dielectric and elas- 
tic properties of the phases a and b with thicknesses x^ and Xb are expressed 
with the high-frequency relative permittivities ea,b and Young's moduli Ej^^, 



respectively. The piezoelectric coefficient for the composite system is then 
calculated as the ratio of the change of surface charge to the force applied 
perpendicular to the surface [ 
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Since there is no intrinsic piezoelectricity assumed in the phases, the coefficient 
dss is written as an effective material property with a superscript 'e'. The term 
on the right-hand side of Eq. (1) containing permittivities is proportional to 
the effective dielectric permittivity of the system, e^ = eaeb(ebXa + eaXb)~^. 
For complex structures, or layered systems with materials having different 
Poisson's ratios, Eq. (1) can not be applied. 

Similarly, applying a numerical model, one should arrive at the same results. 
To this end, we have employed the FE method and performed simulations on 
a double-layer system. In the FE calculations, electrostatic field and stress- 
strain relations are solved simultaneously. Neglecting the polarization of the 
phases, the electric field E(= — V$(r)) distribution is obtained from Poisson's 
equation: 

-V- [5(r)V<l>(r)] = p(r) (2) 

where $ is the electric potential distribution, and e and p are the dielectric per- 
mittivity {e = eeo, Eq is the permittivity of free space and Eq = 8.854 pFm~^) 
and the charge density as space-dependent variables, respectively. For the 
stress-strain relations, Navier's equation is solved. 



-V-T(r) = K(r). 



(3) 
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Fig. 2. Calculated piezoelectric coefficients d^-^ as a function of (a) the relative per- 
mittivity of phase a and (b) Young's modulus of phase a. The symbols (o) and the 

solid lines ( ) represent the results of the FE calculations and the analytical 

expression of Eq. (1), respectively. The dotted ( ) lines represent the permit- 
tivity and Young's modulus of phase b, respectively. The inset in (a) illustrates 
the numerical results for Sa = e^ for various voltages, the slope is the piezoelectric 
coefficient, 246 pmV^^. 



Here T and K are the stress tensor and the body force, respectively. The 
stress T is proportional to the gradient of the unknown displacement u(r) 
and the space-dependent material coefficient c, T = cVu. The coefficient c 
is a function of Young's modulus E^ and Poisson's ratio p of the specific 
medium. The coupling of the electric field to the mechanical stress is achieved 
through the Maxwell stress tensor\l'Ji |20| . 



T'(r) = c(r) Vu(r) - -D(r) • E(r) + E(r)D^(r) (4) 

where T' is the generalized stress tensor. E and D are the electric field and the 
dielectric displacement vectors, and the superscript 'T' denotes a transposed 
matrix. The last term in Eq. (4) is the matrix direct product. Inserting (4) into 
the stress-strain relation of Eq. (3), the general equation for our calculations 
is obtained: 

-V ■ [c(r)Vu(r) - -D(r) ■ E(r) + E(r)D^(r)] = K(r) (5) 

It is worth mentioning that the spatially varying material parameter e and the 
charge distribution p in Eq. (2) are functions of the displacement uJ2l|, |22 . 
However, this coupling is not taken into consideration in this paper. 

For layered structures and for one- dimensional simulations, the non-diagonal 
components in the last term on the left-hand side of Eq. (5) are neglected, and 
the Poisson's ratio v of the materials in question is insignificant. No body-force 
term is considered in the simulations. In the following part, we solve Eqs. (2) 
and (5) with a nonlinear solver based on a commercially available FE software 
package [23^. First two- and three-layer systems with charge distribution(s) 
perpendicular to the apphed field direction (in the y-direction) are consid- 
ered; such a problem is pure one-dimensional. Later the charge distribution 
is considered to be non-uniform and again perpendicular to the applied field 
direction, which makes the problem two-dimensional. We have also assumed 
irregular two-dimensional structures in which the fraction of phase a is de- 
creased in the layer. 



3 Results and discussions 



3.1 Two-layer system 



An application of the numerical model to the geometrical conditions described 
by Kacprzyk et al. (6| confirmed that our approach is valid in one dimension 
with z/a,b ~ and c = E^ . In these simulations, the geometry for the com- 
putation is a layered structure with phases a and b, as shown in Fig. 1. The 



two-layer system is a charged system that produces an external field which 
leads to piezoelectric properties!^. 

The mechanical boundary conditions in the simulations are chosen such that 
the structure is fixed at CD, where it is not allowed to deform in the x- and y- 
directions. At the boundaries AB, AC and BD, the structure is allowed to move 
in the x- and y-directions. The voltage boundary conditions are applied at CD 
and AB with voltages Vcd and Vab, and at AC and BD, symmetry conditions 
are assumed for the static electrical calculations, dV/dx = 0. Moreover, a 
surface charge density of 1 mC m~^ (attainable in experiments |8l. l24| ) is applied 
at the internal EF boundary. The thicknesses of the phases are taken equal 
to each other, Xa = Xb = 1 mm. The layer thickness in the simulations would 
not affect the overall effective properties of the mixture. 

The electric-field-induced mechanical deformation Ay can be expressed as a 
power series. In that case, the electrostrictive Ogg and piezoelectric d^^ co- 
efficients can be estimated from pairs of the position displacements at the 
boundary AB and the voltage difference between the AB and CD boundaries, 
AV = Fab - Vcd. 



= [dl, + al,AV + ...]AV (6) 

where yo is the deformation without the voltage applied and due to the charge 
at the interface EF. We have adapted the first two terms of the series expansion 
at the right-hand-side of Eq. (6) in the numerical data analysis. The numerical 
results are shown as an inset in Fig. 2a. In the inset, the parameters yo, d^^ 
and a^g are 5.65 x IQ-^ [m], 2.46 x 10"^° [pmY^^] and -2.13 x 10~^^[pmV^2]^ 
respectively. Observe that Ogg ■ AV" <^ ^33. The coefficients can be calculated 
with the following relations for the nonlinear, 033 ■ AV ^ d^^, and linear ^33 ^ 
al^- AV regimes, 

al, = AyAV-' [pmY'^] (7) 

dl, = AyAV-' [pmV-i] (8) 

In Fig. 2, the simulation results and their comparison with Eq. (1) are pre- 
sented. The permittivity and Young's modulus of phase b are kept constant, 
eb = 2.5 and E^ = 100 MPa, respectively. First, the permittivity of phase a is 
varied between the free-space value Eq and Aeq while keeping Young's modulus 
of phase a constant, Ej = 1 MPa. Then keeping the permittivity of phase a 
constant, eg = 2.5, Young's modulus of phase a is altered between 100 kPa 
and 1 GPa. These permittivity values are chosen, since most of the polymers 
have relative permittivities between 2 and 12. In addition, traditional ma- 



terials containing voids (1ow-k materials) usually have relative permittivities 
between 1 and 2. Young's moduli of the phases are on the other hand assigned 
on the basis of the effective medium theory [15| for layered structures and the 
measured effective moduli of porous polymers from Neugschwandtner et al. 



24| and of some solid polymers from Mat Web j2 



The results are illustrated in Figs. 2a and 2b. In Fig. 2a, the effective piezo- 
electric coefficient \dl^\ is plotted as a function of the relative permittivity of 

phase a permittivity. The symbols (o) and the solid line ( ) represent the 

results of the FE calculations and the analytical expression Eq. (1), respec- 
tively. There is a maximum of 1^33! when the permittivities of the phases are 
equal to each other, e^ = e^, which is also expected from Eq. (1). This indicates 
that one should tailor or choose materials with matching dielectric properties. 
In Fig. 2b the logarithm of \dl^\ is displayed as a function of Young's modulus 
of phase a for e^ = e^ = 2.5. For a stiff phase a, Ej > E^ , d^^ converges 
to a constant value. As Ej matches the value of E^ , Ej ^ E^ there is no 
piezoelectric activity, log(ii^J) = 8. As the stiffness of phase a becomes lower 
than that of phase b, \dl^\ increases. For Ej <^ E^ there is a linear relation 
between Ej and \dl^\ in the log-log graphs. The simulation indicate that there 
are are no nonlinear effects of the applied voltage on the deformation in these 
two-layer simulations, so that 0^3 ?« in Eq. (6), see inset in Fig. 2a. All 
these observation are also expected from Eq. (1), which is also presented in 
the figures with solid lines ( ). 

Next, we examine the influence of the charge distribution at the interface on 
d^g. To this end we introduce a volume charge distribution p(x, y) which is only 
varied in the ^/-direction. The validity of Eq. (1) for such cases is discussed 
below. The charge distribution is generated with a delta sequence J26|, 

My) = nn''[i + nW (9) 

where 

\im J Uy)dy = I. 

In Eq. (9), n is a shape parameter. For charge distributions in the form of 

p{x,y) = po(j)n{y); 

the volume charge distribution becomes surface charge distribution p = po6{y) 
as n -^ 00. 

In Fig. 3a the calculated \dl^\ is presented as normalized with the volume 
charge density po/Qn for different shape parameters n, where 

Qn= / Po • (t>n{y)dy. 





260 


1 — 1 




1 


?4n 


> 




b 


??n 


Q. 






200 


O 




"^^o 


180 


Q. 




X 

CO 


160 


CO 




^ 


140 




T 










CM 




1 

E 


2 


O 


1 


b 




n 


a- 








D 


-1 


o 






_*«# 


V (b) ■ 


-^--^^yy 


\5n^-^ 


— = 5^^^/yy 


\\ \ ^^ -^— — = — 




\ ~\ ^^^ ^^^^ ' — .^_ - 


^ — " ^^^ ^^ ^ 


\ ^^ ^^^ ^ — — ^ 


^=^^^ 


^^^=^ 



-100 



-50 





y[^m] 



50 



100 



Fig. 3. (a) Calculated absolute values of piezoelectric coefficient {d^^] normalized 
to the actual delta surface charge density po and the total assumed volume charge 
density Qn = j p{y)dy. The symbols (+) are results from the FE calculations, the 

solid line ( ) represents the actual value from Eq. (1), \d^^\ = 24:7. b pmV^^ for 

a charge sheet, (b) Charge density distributions at the interface for 3 < n < 6 in 
Eq. (9. 
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Fig. 4. Computation domain for the three-layer system. 



The total volume charge density is not equal to po because of the finite n. 
However, as n gets closer to 10^, Eq. (9) can be used to represent a charge 
sheet. In Fig. 3b volume charge distributions are illustrated for the values of n 
used to calculate cigg in Fig. 3a. As the charge penetrates inside the material, 
forming a true space charge, the piezoelectric response of the composite system 
becomes weaker compared to a uniform sheet-like surface charge distribution. 



Table 1 

Minimum and maximum values of the parameters used in Eqs (6) for various sim- 
ulations. Order means the layer sequence, soft-rigid-soft (SRS) and rigid-soft-rigid 
(RSR). The parameters varied in the simulations are also presented. 



Order/ 
Parameter 


a|3 [fmV-2] d|3 [mV-^] 
min max min max 


2/0 [M-m] 
min max 




Two-layer system 


/fa 

—lEl 


204p 250p 
-2.00p 2.49n 


4.72 5.70 
0.07 56.5 


bipolar charge system 



RSR/ea 

SRS/ea 
RSR/S3Y 

SRS/£'a 



RSR/ea 
SRS/ea 

RSR/^a^ 

SRS/^a^ 



2.27 2.79 -13. 9f -11.4f 

2.27 2.80 -245p -200p 

0.03 27.7 -2.47n 2.20p 

0.03 27.7 -2.20p 2.47n 



mono-polar charge system 



2.27 2.79 11.3f 14.0f 

2.27 2.79 -13. 9f -11.4f 

0.03 27.5 -0.13f 0.15p 

0.03 27.7 0.14p 0.13f 



0.23 


0.25 


4.66 


5.61 


0.06 


55.5 


0.06 


55.5 



13.98 


55.9 


0.23 


0.25 


0.02 


223.7 


0.02 


0.32 



3.2 Three-layer systems 



In order to enhance the efficiency of an electro- mechanical transducer, it is 
possible to stack layers of different materials. Such a system can be the three- 
layer binary structure ('a sandwich') which allows to employ and investigate 
different charge polarities. We therefore apply the numerical method to a 
three-layer binary system as shown in Fig. 4. In these simulations, (i) bipolar 
and (ii) mono-polar charge systems are assumed. In a bipolar system, the 
charge polarities at interfaces EF and GH in Fig. 4 are opposite to each other. 
In a mono-polar system, on the other hand, they have the same polarity. (The 
main difference between bi- and mono-polar charge systems is that bipolar 
charge systems do not contain unbalanced internal charge). The deformations 
Ay are analyzed as functions of the applied voltage difference Vab — Vcd 
with Eq. (6). The deformation Ay of a bipolar charge system is dominated 
by the contribution of the linear term, ligg. The deformation of a mono-polar 
charge system, on the other hand, shows negligible linear dependence and 
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Fig. 5. Thickness change y — yo in a three-layer structure for (a) bi- and (b) 
mono-polar charge distributions as a function of the voltage difference between 
the electrodes Ay = Vab — Vcd- The symbols are the values calculated with the 

FE method. The solid lines ( ) are the best linear and quadratic fits. In (b) the 

fitted curve is separated into its linear ( ) and quadratic ( ) parts, in order 

to show the electrostrictive and piezoelectric contributions. 

the quadratic term, Ogg is important indicating 'electrostrictive phenomena'. 
The calculated deformations for these two systems are shown in Fig. 5. The 
deformation with the same applied voltage is 1000 times higher for bipolar 
than for mono-polar charge. 

In the simulations, the surface charge densities are ±1 mC m^^ and 1 mC m^^ 
for bi- and mono-polar charge distributions. Again, similar to the two-layer 
system, we first alter the permittivity of the phase a keeping Young's moduli 
of both phases and the permittivity of the phase b constant. Later, Young's 
modulus of phase a is varied while all other material parameters are kept 
constant. The total thicknesses of the phases are equal to each other x^ = 
x\j. There are two different combinations of interest: (i) phase a in Fig. 4 is 
rigid and the structures is soft-rigid-soft (SRS) or (ii) phase a is soft which 
leads to a rigid-soft-rigid (RSR) structure. The mono-polar-charge cases do not 
produce significant electro-mechanical activity comparable to the bipolar ones. 
Experimental observations performed on some prepared RSR structures, have 
verified the numerical simulations such that the RSR structures do not yield 
any observable piezoelectric activity. Moreover, mono-polar-charge samples do 
not show any measurable electro-mechanical activity. 

The results of the simulations are illustrated in Fig. 6. In the figures, the 
values are normalized respective to a maxima of the absolute parameter val- 
ues, maxlalgl, maxidggi and max||/o|- The normalization values are listed in 
Table 1 as min and max. In the table, the rigid and soft phase order of the 
structures are presented together with the parameter that is varied. While 
the total thickness change behaves differently, the normalized electrostrictive 
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Fig. 6. Calculated (a) electrostriction a^^ (b) piezoelectric coefficients dig and (c) 
deformation yo due to interface charge for altered permittivity of phase a. (d) Elec- 
trostriction a|3 (e) piezoelectric coefficients d^^, and (f ) deformation yo due to inter- 
face charge for altered Young's modulus of phase a. The parameters are normalized 
with respect to their max or min values, which are presented in Table 1. The filled 
and empty symbols represent the mono- and bipolar charge cases for the simula- 
tions. The solid lines ( ) in (b) and (e) are obtained from Eq. (1) with interface 

charge p = 970 and 992 [xCm"^, respectively. 

and piezoelectric coefficients had the same shape as functions of the relative 
permittivity of phase a. 

In Figs. 6b and 6e, solid lines ( ) represent the curve-fitting results obtained 

with Eq. (1) in which the surface charge p is used as a free parameter. It is 
found that ±1 mC m~^ on the interfaces of a three-layer system yields the same 
results as a two-layer system with approx. 1 mC m~^ at the interface. This is 
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Fig. 7. (a) Calculated piezoelectric coefficient |d|3| (•) as a function of difference in 

total charges at boundaries EF and GH. The solid line ( ) is Eq. (1) with half 

the apparent charge, (b) Relative error in percent. 



due to the symmetry plane between the two charge layers that can be assumed 
as in the method of images 
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Since there is no external field in a three- 
layer system with bipolar charge-it is neutral-it is to be preferred for non- 
electrostatic industrial applications. Except for the results of the RSR/eg case, 
the other three bipolar charge systems result in high piezoelectric coefficients, 
^33 > 200 pmV~^. It is interesting that the max|(i33| and min|(i33| values 
obtained for the RSR/ea and SRS/ea cases are not comparable. The RSR 
structure is not suitable for any electro-mechanical application. The reason 
for the very low piezoelectric coefficients in the RSR case is that the structure 
expands or contracts in response to the internal charge distribution, even in 
the absence of an applied voltage-|/o in Table 1. After this deformation there 
is no (or very little) space for the structure to respond to the applied voltage. 
It is observed that a ^33 value on the order of nmV~^ can be obtained if 
materials with low Young's modulus {E^ < 1 MPa) are employed such as 
non-conducting polymer foams and natural corkJ2 



As mentioned previously, the mono-polar charge systems exhibit very little 
electro-mechanical activity when the charges at the surfaces are equal. How- 
ever, at high applied voltages, the electrostrictive effect is significant. To have 
the same magnitude of the piezoelectric coefficient in a two-layer system one 
should introduce a surface charge of about InC m~^ at the interface, which is 
extremely low. If the surface charges are not equal, but have the same sign, 
the response of the structure can be tailored to be linear in the applied field 
for specially selected surface charge values. This is indicated in Fig. 7. When 
the unbalanced charge of the structure is close to zero, (pef — Pgh)/2 ~ 0, 
the analytical solution and the numerical results deviate, due to the numerical 
error. This is an interesting case, since in phases with significant conductivity. 
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charges on the surfaces would recombine. This process would directly affect 
the electro-mechanical activity and lower d\ 



e 
33- 



3.3 Two-dimensional system 



3.3.1 Non-uniform charge density at the interface 

As an illustration for a two-dimensional solution, we again consider a two-layer 
system as in Fig. 1. However, this time we just alter the charge distribution 
at the interface EF in the x-direction. 

p{x) = pomr-^[l + r?x\^ (10) 

where n defines the shape of the distribution — for n — > surface and n ^ oo 
line charge distributions are obtained — , and po is the charge amplitude. The 
material properties are taken as Ca = eb = 2.5 and E^ = 1 MPa and E^ = 
100 MPa. The boundary conditions are assigned as in the previous problem 
(see § 3.1). Since the problem is now in two-dimensions the Poisson's ratio of 
phases u^ and u^ influence the resulting effective piezoelectric coefficient d^^. 
In Fig. 8a, the results obtained for d^^ are illustrated as a function of n for 
various Poisson's ratios. For small values of n (n < 10^) and Ua = u^ <^ 0.1, rfgg 
is the same as the one-dimensional simulations, a uniform charge distribution. 
However, as the charge distribution is altered by changing n, the resulting ^33 
increases for n > 10^ and approaches a constant for n > 10*^. These higher 
values of ^33 are expected due to the electric field distribution in the system, 
which is enhanced at the tip of the charge distribution (sheet-like discontinuous 
surface charge in the x-direction). The enhancement is higher as the charge 
distribution approaches a line charge distribution, the piezoelectric coefficient 
being approx. 10% higher than that of uniform sheet charge distribution. In 
Fig. 8a and 8b show that ^33 decreases with increasing Poisson's ratio, when 
the latter approaches 0.5. 

Finally, as an illustrative example in Fig. 9a the stress distribution and elec- 
trical potential are shown. The voltage difference between the AB and CD 
boundaries is V. The Poisson's ratio of the phases are as follows, u^ = 0.33 
Uh <^ 0.1. In Fig. 9b the charge distribution at the interface and the defor- 
mation on the boundary AB are presented. The deformation is purely due to 
the interfacial charge, and it is larger than the previous cases considered. The 
deformation is localized as the charge distribution. 
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Fig. 8. (a) Calculated absolute values of piezoelectric coefficient \d^^\ normalized to 
the actual two dimensional charge density p{x) for various Poisson's ratio of phases 
(u = u^ = z/b); (□) ly < 0.1, (x) u = 0.1, (+) u = 0.33, (*) u = 0.5 and (A) 

fa *C 0.33 i^t, = 0.1. The solid line ( ) is the value calculated in § 3.1 for uniform 

surface charge distribution, (b) Influence of Poisson's ratio on dig for u^, = Ub <^ 0.1. 
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Fig. 9. (a) Mechanical stress (von Mises stress) in gray-scale and voltage distribution 
in contour plot (for AV = V). The structure is deformed, the maximum and 
minimum deformations are marked with symbols ( x ) , and the values are normalized 
to thickness of the whole structure, AC. (b) Charge density at the interface EF and 
deformation of the boundary AB. The charge density is calculated for n = 10 and 
Po = 1 mCm~^ in Eq. (10). 

3.3.2 Irregularly shape sample 

In this section we have considered an irregular shape as illustrated in Fig. 10a. 
Phase a is considered to be changing its width, |AP|. In the simulation, the 
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Fig. 10. (a) The geometry used in the calculations for the irregular shape. Phase 
b is considered to be soft, the length PB is used as a parameter in the calcula- 
tions, X = |PB|/|CD|. The dotted line between EF show the charge interface, (b) 
Calculated absolute values of piezoelectric coefficient 1^33! for various PB ratio and 
Poisson' ratio when the charge distribution of 1 mCm^^ is considered in the EF 
interface; (■) i^^ = 0.33 and i^b < 0.1, (+) fg < 0.1 and i^^ < 0.1, (•) i^^, < 0.1 
and Ub = 0.33, (v) Va = 0.33 and vt, = 0.33, (o) u^ = 0.5 and v^ = 0.33. When only 
charge is considered in the EO the piezoelectric coefficient is not altered significantly, 
(A) Ua = 0.5 and i^^ = 0.33. 



0.002 




Fig. 11. Mechanical stress (von Mises stress) in gray-scale and voltage distribution 
in contour plot (for AV = V). The structure is deformed, the maximum and 
minimum deformations are marked with symbols ( x ) , and the values are normalized 
to the thickness of the whole structure, AC. The Poisson's ratio of phase a is 0.33. 



same consideration as in the previous cases are considered, €3 = e^ = 2.5 and 
EJ = 1 MPa and E^ = 100 MPa. The charge is assumed to be on the interface 
EF with the charge density 1 mCm~^. The results are plotted in Fig. 10b for 
various Poisson's ratio of phases and fraction of phase a, 1 — a; where x is the 
length of PB. It is clear that the Poisson's ratio of the soft phase (phase b) 
influences the piezoelectric coefficient \dl^\. Moreover, as the fraction of the 
rigid phase is lowered, a; — > 0, the piezoelectric activity is decreasing linearly 
with the fraction of phase b, which is expected. In addition the decrease in ^^3 
with respect to x is proportional to the Poisson's ratio of phase b. In Fig. 11 
the stress and the electric potential distributions are presented as gray-scale 
and contour plots. 
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4 Conclusions 



Numerical simulations on the electro-mechanical properties of layered struc- 
tures containing charges are reported. The simulations take the coupling of 
the electrical and mechanical stresses into account by means of the Maxwell 
stress tensor. The results are compared with an analytical model, and it is ob- 
served that there is good agreement for two-layer systems. The piezoelectricity 
in three-layer systems shows a strong dependence on the polarity of interface 
charges. It is concluded that bipolar charge systems clearly yield higher piezo- 
electric coefficient values than mono-polar charge systems. The mono-polar 
charge cases are dominated by 'electrostrictive effects'. 

In this paper, we have presented a way of calculating electro-mechanical ac- 
tivity of layered electrets with uniform and non-uniform charge distributions. 
The numerical model is also applied to a two-layer system with discontinuous 
interface charge layer. In such a simulation unlike the one-dimensional cal- 
culations the Poisson's ratio of the phases become significant. High Poisson's 
ratio lead to low the electro- mechanical coefficient. As the sheet of charge at 
the interface is modified toward a line charge distribution the piezoelectric co- 
efficient is increased approx. 10%. The two-dimensional structures considered 
have also illustrated the stress distribution in the constituents are of impor- 
tance, the stiff phase is mechanically stressed due to electromechanical activity 
in the soft phase. We might conclude that the presence of excess charges (space 
charge) in heterogeneous systems is crucial for aging of the materials. 

As a concluding remark, it is illustrated that in heterogeneous materials, fluc- 
tuations in the local electric field could lead to mechanical deformations even 
if there is no intrinsic piezoelectricity in the phases. The considered cases show 
that the deformation follows the actual charge distribution at the interface, 
and distribution of charges plays an important role in materials with irregu- 
larity. 
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